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Abstract 

We report results on the nucleon structure obtained from the lattice quantum chro- 
modynamics calculation. These include the axial, electromagnetic, ttNN, and scalar 
form factors. The calculation is carried out at /3 = 6 on a 16'^ x 24 lattice with 24 
quenched gauge configurations. The chiral limit results are extrapolated from several 
light quark cases. For the disconnected insertion (sea-quark contribution), we used 
the stochastic estimation with the Z2 noise to calculate the diagonal and off-diagonal 
traces of the inverse matrices with a size of 10^ x 10^. It is found that the Z2 noise is 
the optimal choice and its comparison with the Gaussian noise for our quark matrix is 
given. For the sea-quark contribution, we report results on the strange condensate in 
the nucleon and the nNa term. 



1 Introduction 

The Understanding cannot See. 
The Senses cannot Think. 

By their union only can Knowledge be produced. 

— Immanuel Kant 

Imagination is more important than Knowledge. 

— Albert Einstein 



The above quotes can be used to describe the two well recognized approaches to scien- 
tific research — Experimental and Theoretical. Particularly true to the quotes is the field 
of strong interaction where there has been an intense interplay between theory and exper- 
iment. The advent of quantum chromodynamics(QCD), a non-abelian gauge relativistic 
quantum field theory of quarks and gluons (the constituents of protons and neutrons), in 
the early seventies has offered the most promising description of strong interaction dynam- 
ics and internal structure of hadrons both in logical consistency and in scope It is widely 
accepted that it is the correct theory of strong interaction. This acceptance is based on a 
substantial amount of experimental data accumulated and refined over the last two decades. 
It is particularly so for the inclusive scattering processes at large momentum transfer where 
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the perturbative QCD is applicable. On the other hand, confirmation of QCD as the funda- 
mental theory of strong interaction has been somewhat hampered by the lack of analytical 
calculations of the low frequency modes due to the inherent non-perturbative nature of the 
theory. 

The invention of lattice-regularized QCD||2[ with Monte Carlo methods]^ holds the 
promise of circumventing the difficulty of obtaining analytic results via numerical simula- 
tion. In the lattice-regularized gauge theory, the continuous space-time is discretized by a 
finite lattice with a periodic boundary condition. As a consequence, the infinite dynamic 
degrees of freedom associated with the space-time is reduced to a finite number which makes 
numerical analysis feasible. The recent advancement in supercomputer technology and the 
DOE Grand Challenge Award have made it possible to carry out large scale Monte Carlo 
simulations to calculate the hadronic structure and interaction in the framework of lattice 
QCD. This allows us to compare the lattice results with the experimental data in a quali- 
tative and semi-quantitative manner which is very valuable in understanding physics from 
first principles free from the uncertainty of models. 

In the present paper, we shall report results on the study of the nucleon structure which 
includes the electromagnetic, axial, pseudoscalar, and scalar form factors and compare 
them with the known experiments. From the pseudoscalar form factor, we extract the pion- 
nucleon-nucleon form factor. We also report results on the sea-quark contributions with 
disconnected current insertions. This entails quark-loop calculations involving diagonal 
and off-diagonal traces of the inverse quark matrices with a dimension 10^ x 10^. This is 
carried out with the stochastic estimation algorithm with the noise. We show that the 
Zi noise is actually the optimal choice and a direct comparison with the Gaussian noise is 
given. From this calculation, we obtain the strange quark condensate in the nucleon. 

As we shall see the above computer simulation has already been very helpful in bridging 
between the theory and experiment via ah initio calculations and is doing a job better than 
the models, this third branch of scientific research will only get more mature when both the 
computer hardware and algorithm improve with time. To put the numerical simulation in 
perspective, we maintain: 

Imagination may not be relevant. 
Knowledge cannot predict. 

By simulating Imagination and testing against Knowledge 
only can Reality be re-created. 

2 Formalism and Numerical Steps 

The Euclidean path integral formulation of quantum field theory is employed to calculate 
the Green's function. In this way, the quantum field theory becomes a problem in classical 
statistical mechanics. Physical observables are then extracted from the statistical correlation 
functions of composite operators built on the fundamental dynamical variables in quark and 
gluon fields. In the following, we present a synopsis of the formalism and related numerical 
steps. 
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The correlation functions have the following generic form 



(1) 

where Z is the partition function. Sq is the action of the gauge link variable U and Sp 
the fermion action which is bilinear in the quark field variable ^ and can be written in the 
form Sp = '4)M\U]iIj. Since both the fermion action Sp and the operator product O1O2... 
are bilinear in tp and ^, the fermion integral over the anti-commuting Grassmann variables 
■0 and ■0 can be done analytically to give 

< Oi{U,i),^)02{U,i),^)... >= ^ J [dUje-'^'^ det M[U]Tr[M-^[U]M-^[U]M-^[U]...T]. 

(2) 

The path integral in eq.(§) is an ordinary multiple integral over the group manifold, in 
this case SU (3) group. One might think of approximating the integral by summing over 
a sufficiently dense set of mesh points, say 4 per variable. However , this "modest" goal 
is well beyond the capacity of the existing supercomputers: for a (small) 10^ lattice, it 
would require a calculation on 4^20,000 _ -|^q192,659 points. Instead, one may use the im- 
portant sampling technique to select a comparatively small subset of gauge configurations 
{Ui,U2, ■■■,Uiy) among the "important" ones, such that the probability of occurrence of a 
given configuration Ui in this ensemble approaches the desired distribution 

[dU]e-^'^ det M[U], (3) 

for — > 00. Then the quantum average are given(for sufficiently large N ) by averages 
taken over the sample: that is 

1 ^ 1 
< O1O2... >= -Y.Tr[M-\U,]M-\U,]...T]±0{-^). (4) 



There are three numerical steps to lattice gauge calculations. The first step is to prepare 
an ensemble of gluon configurations {Ui} in the vacuum. This is done using the Cabibbo- 
Marinari[P pseudo-heatbath algorithm, a Monte Carlo method which generates sets of link 
matrices {Ui} according to the distribution in eq.(^). In principle, the probability distri- 
bution includes a factor of det(M[?7]) which represents the eff^ect of quark-antiquark polar- 
ization in the vacuum. In the present simulations, we omit this factor as an approximation 
since it would require more computation than we can do. We will discuss the validity of this 
approximation when we come to the results of the calculations. Omitting this determinant 
is referred to as the quenched approximation. 

The second step is to calculate the quark propagator matrix M^^[U] in the prepared 
gluon ([/) background. This is done with the conjugate- gradient algorithm[p. This is an 
iterative scheme which is capable of giving the desired solution to arbitrary accuracy. We 
typically find a few hundred iterations to be sufficient. 

The last step is to assemble the quark propagator matrices M~^[?7j] and the the matrix 
T[Ui] defined in eqs. (|l|) and @ to evaluate the correlation function in eq.(^) from which 
we extract the physical quantities with statistical analysis. 
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3 Form Factors of Nucleon 



3.1 Formalism 

Electron and neutrino scatterings off the nucleon target are useful probes of the nucleon 
internal structure by measuring its electromagnetic and axial form factors. These form 
factors can be simulated in the lattice gauge calculation and compared with experiments 
directly and thereby serve as a test to the Monte Carlo calculation and the theory as well. 
The irNN coupling form factor is a fundamental quantity in the pion-nucleon and nucleon- 
nucleon dynamics. Yet it is not experimentally accessible directly and is poorly known 
theoretically. From the calculation of the pseudoscalar form factor, we extract the ttNN 
form factor with pion pole dominance. In this case, it is a prediction from the fundamental 
theory and will have great help in settling the large uncertainty over the shape and the 
cutoff of the ttNN form factor. Furthermore, there have been sprouting interest in the 
sea-quark contributions to the nucleon structure, for example, the irNa term, the strange 
condensate in the nucleon and the flavor-singlet axial charge g\. The recent deep inelastic 
scattering experiment obtained a very small g\ which suggests that most of the spin of the 
proton does not come from the quark spin which is a great surprise and has been dubbed 
the "proton spin crisis" . Since no model we know has a handle on this problem, we expect 
the answer to come from the lattice QCD calculation. 

The nucleon can be obtained from the following two- and three-point correlation func- 
tions |, g § 

GNN{t,p) = Y.e-'P-^ < 0\N^{x)No,mO > . (5) 

X 

G7Vj7v(ti,t2,P,q) = E e^^-"^-^P-"-r"^ < 0|iV^(x2)i;.(xi)iV,(0)|0 >, (6) 

Xl,X2 

where represents the proton or neutron interpolating field with the Dirac component a, 
r is a 4 X 4 matrix in the Dirac space and represents various bilinear quark currents in the 
form of '^S^'. By inserting complete sets of states and setting t2 >> ti >> a, the lattice 
spacing, we obtain the electromagnetic form factors GE{q'^) and GMi'f'): axail form factor 
9a{(1^)-: pseudoscalar form factor gp{q'^), and scalar form factor gs{(f ) from the ratios of the 
three point functions involving corresponding currents to two point functions. The results 
corresponding to the physical quark masses are obtained from the extrapolation from those 
with heavier quarks to the chiral limit where the pion mass is zero. 

3.2 Numerical details 

Our results are based on 24 quenched (ignoring the determinant in eq. (^) gauge configu- 
rations on a 16^ X 24 lattice and were calculated using the Monte Carlo Cabibbo-Marinari 
pseudo-heatbath algorithm[p. The 5C/(3) fundamental Wilson action was used with peri- 
odic boundary conditions and the coupling constant was set at [i = 6.0. The gauge field was 
thermalized for 5000 sweeps from a cold start and 24 configurations separated by at least 
1000 sweeps were saved. For the quarks we use periodic boundary conditions in the spatial 
directions and "fixed" time boundary conditions, which corresponds to setting the quark 
couplings across the time edge to zero. The origin of all quark propagators was chosen to 
be at lattice time site 5; the secondary zero momentum nucleon source was fixed at time 
site 20. We expect that these positions are sufficiently far from the lattice time boundaries 
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to avoid nonvacuum contaminations. We used the red-black pre-conditioned conjugate- 
gradient algorithm for the quark propagator. For our convergence criterion we demanded 
that the absolute sum of the squares of the quark propagators be less than 5 x 10~^ over 5 
iterations. As a check of the nucleon secondary source, we verified current conservation for 
t2 > ti > to 0(10-4). 

3.3 Stochastic Estimation with Zi Noise 

Our present space-time lattice with merely the size of 16'^ x 24 gives a quark matrix of 
the dimension 10^ x 10® including the spin and color degrees of freedom. While it is 
durable to calculate the quark propagator, i.e. M-^(x,0) for a point source S at with 
a reasonably small quark mass (e.g. a fraction of the strange quark mass) on today's 
supercomputers, the quark propagator M-^(x,y) from any point to any point is certainly 
unattainable. For calculations of the 2-point functions (eq. (^) and 3-point functions (eq. 
(^ with connected insertions for flavor non-singlet currents, one can get by with the help 
of translational symmetry and uses only M-^(x,0). But there are cases where one can not 
rely on such a help. These include the calculations of quark loops which are space-time or 
space integrations of the fermion propagators. Examples of interest in QCD include the 
quark condensate and the topological susceptibility with the fermion method, flavor-singlet 
meson masses which involve disconnected quark loops in the two-point functions, notably 
the U(l) problem., and the vrA^ a term and the proton spin problem which involve quark 
loop contributions in the three-point functions. 

Instead of waiting for the advent of more powerful hardware, we have employed the 
stochastic algorithm with the Z2 noise to estimate the quark loops |p . Stochastic approach 
to estimating the inverse of an A'^ x matrix M entails the introduction of an ensemble of 
L column vectors 77 = 77^, ...,?7^ (each of dimension A x 1) with the properties of a white 
noise, i.e. (r/j) = 0, = The expectation value of the matrix element can be 

obtained by solving for Aj in the matrix equations MX = rj with the L noise vectors r] and 
then take the ensemble average with the j-th entry of r] 

E[Mr^^] = {n,X,) = J2 Mr^Hii.Vk) = M-.\ (7) 

k 

which is the matrix element itself. In fact, it has been shown recently ||l^ that the 

variance of a inverted matrix element due to the stochastic estimation is composed of two 
parts 

^ Mi 

Whereas the second part is independent of the kind of noise used, the first part is propor- 
tional to the square of the diagonal error C2 = \J jj J^iiiViVi) ~ 1)^- Since Z2, or Zn for 
that matter, has no diagonal error, i.e. C2 = 0, it produces a minimum variance. Other 
noises will have larger variances due to the non- vanishing C2. For example, C2 = \/2/L for 
the Gaussian noise with large and independent configurations. 

The question remains as to whether the Z2 noise is superior than the other noises on a 
practical footing with reasonable small L. To answer this question, we have tested the use- 
fulness of the Z2 noise against the Gaussian noise for small noise configurations L (L < 100) 
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and smaller quark masses. Plotted in the left column in Fig. 1 are the accumulated averages 
for the estimate of the real part of the diagonal trace -Rerr^^M~^(x,x)/y as a function 
of L for the and Gaussian noises with Wilson hopping parameters k = 0.148,0.152 and 
0.154 for a gauge configuration. With Kc = 0.1568, k = 0.152 corresponds to the strange 
quark mass and k = 0.154 corresponds to a mass about half of that. The straight line is the 
result from the stochastic estimate at L = 300. We see that the estimate from the Z2 noise 
approaches the value at L= 300 (assumed to be the asymptotic value) faster than that of 
the Gaussian noise. The corresponding jackknife errors for different L (the right column in 
Fig. 1) show that the Z2 noise is consistently better than the Gaussian noise by a factor of 
two for all three k's. This means that to achieve the same level of accuracy for the diagonal 
trace, one needs statistics for the Gaussian noise about 4 times as much as for the Z2 noise. 
We have also examined the near diagonal trace related to the point-split axial current. The 
result is similar to that of the diagonal trace @. 

We have employed the Z2 noise to calculate the quark loops in the presence of the 
nucleon with L = 200 to 300 only. Some of the results will be reported in the next section. 
Compared with the brute force approach of inverting the whole quark matrix, we have saved 
the computer time by as much as a factor of 6000! 

It is worthwhile noting that the stochastic estimation is particularly successful for the 
trace (denotes as Re^^ in Fig. 1). This is due to the translational, color and spin sym- 
metries. As a result, the error is proportional to l/\/iV where N is the dimension of the 
matrix. With N = 1.18 x 10^ in our case and the ratio Efe^ii^fcl^lVl^f/]^ = 0-8, we 
predict the error to signal ratio to be \/1.5 x 10~^ from eq. for L = 1. This agrees well 
with the numerical calculation shown in Fig. 1. Given this level of accuracy, it is feasible 
to apply the stochastic method to the calculation of the determinant, the eigenvalues, and 
the eigenvectors of the matrix M which might not be feasible with other algorithms. 



4 Results 

We present our results in three groups. The first group includes the electromagnetic and 
axial form factors where the experimental results are available. We shall compare them 
with our calculation directly. The second group consists of the ttNN form factor deduced 
from the pseudoscalar form factor where is no direct experimental result available. The 
third group involves flavor-singlet quantities which require the quark loop calculation with 
disconnected insertions like the nNa and the strange condensate in the nucleon. 



4.1 Electromagnetic and axial form factors 

In Fig. 2, we plot the isovector axial form factor extrapolated to the chiral limit for the 
local current (L.C.) and the point-split current (P-S.C.) in comparison with the exper- 
imental result. In doing so, we have used the calculated nucleon mass to set the scale 
for the momentum [Q]. In addition, we also show the calculated electric form factor Ge 
extrapolated to the chiral limit and the corresponding experimental results. The experi- 
mental QAil^) has been measured in neutrino-neutron scattering and pion electroproduc- 
tion. The neutrino data gives a good fit in the dipole form up to = 3GeV^/c^ |11|, i.e. 



9A{q ) = 9A(0)/(l-g /M^)^, with the axial vector coupling constant ^^(0) = 1.254±0.006 
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and Ma = 1.032 it 0.036GeV (world average). The new data from Brookhaven E734 exper- 
iment [12| gives a value Ma = 1.09 it 0.03 it 0.02 which is higher than the previous world 
average. Our fit of the axial form factor to the dipole form yields Ma = 1.09 it 0.05GeV 
for the P-S.C. and Ma = 1.06 it 0.04GeV for the L.C. which are closer to the new 
experimental fit of the Brookhaven E734 data [|l^]. Similarly, the fitted dipole mass 
for Ge is Me = 0.857 ib 0.036GeV which is close to the experimental dipole mass of 
0.828 lb 0.006GeV |l^. Comparing to the experimental value qa = 1.254(6), we find that 
the calculated qa = 1.20(11) and 1.18(11) from the P-S.C. and L.C. are 4% and 6% smaller 
respectively. 



4.2 ttNN form factor and the induced pseudoscalar form factor 

The irNN form factor g-wNNil^) is a fundamental quantity in the low-energy pion-nucleon 
and nucleon-nucleon dynamics. Many dynamical issues like the vrA^ elastic and inelastic 
scattering, NN potential, three -body force (triton and He^ binding energies), pion photo- 
production and electroproduction all depend on it. Similarly, the pseudoscalar form factor 
is important to testing low-energy theorems, chiral Ward identity and the understanding 
of the explicit breaking of the chiral symmetry. Yet, compared with the electromagnetic 
form factors and the isovector axial form factor of the nucleon shown in the last section, 
the pseudoscalar form factor gp{q^) and the ttNN form factor g-wNNio^) are poorly known 
either experimentally or theoretically. Notwithstanding decades of interest and numerous 
work, the shape and slope of ^TrAfAf (9^ ) remain illusive and unsettling. Upon parametrizing 
g-nNNio^) in the monopole form or the dipole form, the cutoff mass can differ as much as 
a factor of 3 in different models. 

The ttNN form factor (shown in Fig. 3) is obtained from the pseudoscalar form factor 
with pion pole dominance [^J]. We fitted it with both a monopole form and a dipole form. 
We found that the monopole form with At^nn = 0.75±0.14GeV agrees with the Goldberger- 
Treiman relation at = and the dipole form does not. Extrapolation of the monopole 
fit to q'^ = m\ gives g-wNN = 12.7 it 2.4 which agrees well with the phenomenological value 
of 13.40 ± 0.17. g-nNN{'^) = 12.2 ± 2.3 which agrees with the GT relation. 

This is a good example where the important physical quantity is not directly accessible 
experimentally or reliably obtainable from the models, lattice calculation can fill the void 
and give a reliable prediction. 



4.3 ss in the nucleon and nNa term 

The matrix element (A^|ss|A^) and the i^Na term requires the quark loop calculation with 
disconnected insertions. 

For the disconnected insertion of the current J{x,t) = ^{x,t)T^{x,t), we calculate 
the ratio of the three- and two-point functions. 

^(^^^tio)) -E('^(^>^)) -^t»a const + t{N\^Tm)ci^s. (9) 

Here is the nucleon interpolation field and x is summed over to obtain the forward 
matrix element. Hence, the matrix element can be obtained as the slope from the above 
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ratio. Since we use the fixed boundary condition for the quark field in the time direction, 
r is summed to 4 steps away from the time-boundary in both ends to avoid the boundary 
effect. 

Plotted in Fig. 4 are the ratios defined in eq. (^) for quark masses ranging from the 
charm (Wilson k = 0.120) to strange (k = 0.154). We see clearly that the slope becomes 
larger when the quark mass becomes smaller. Extrapolating to the chiral limit, we obtain 
{N\uu + dd\N)(iis = 5.36 it 1.04. This is very large, about twice of the connected insertion 
{N\uu + dd\N)con = 2.86(6). If we take the quark mass to be 5 MeV, the calculated vrA^fi 
term is then 41.1 it 5.2 MeV. This is quite in agreement with the 45 MeV obtained from 
vrA^ scattering and is also consistent with a similar calculation on a 12^ x 20 lattice with 
the volume source [15|. 

We have also calculated (A^|ss|iV). In this case, the quark mass in the loop due to the 
current self-contraction corresponds to the strange with k = 0.154. The valence quarks in 
the nucleon interpolating filed are then extrapolated to the chiral limit. For the strange 
condensate, we obtain (A^|s,s|A^) = 1.72 it 0.28. If the strange quark mass is taken to be 
130 MeV, then the strange will contribute ~ 224 it 36 MeV to the nucleon mass. This 
constitutes an appreciable percentage of the total nucleon mass. 



5 Computational Aspects and Conclusion 

We have demonstrated in this paper that it is feasible to perform ab initio calculations of the 
hadronic structure directly from the fundamental field theory — quantum chromodynamics 
without having to rely on models. Notwithstanding the fact that the present lattice size 
is modest ( e.g. 16^ x 24), the quark masses used are heavy compared to the physical 
situation, the quenched approximation is used, and the Z2 noise algorithm saved us 6000 
times of the computer time for the quark loop calculation, present calculations still took 
substantial supercomputer resources to perform. For example, the project on the quark 
loop calculation took 2000 C-90 CPU hours at 1 GFLOP/processor speed to accomplish. 
For the 16^ x 24 lattice, the memory requirement is 20 Megawords. 

In order to test QCD as the fundamental theory and make predictions to compare with 
experiments beyond doubt, we need to push the calculation to a much larger lattice, a 
smaller quark mass and with the full description of the dynamic fermions (non-quenched 
approximation). Each of these directions will demand one or more orders of magnitude in 
speed and memory than the present calculation requires. 

This work is partially supported by DOE Grant DE-FG05-84ER40154. The author 
is indebted to S.J. Dong, T. Draper, W. Wilcox, and C.M. Wu who collaborated on the 
projects presented here. He is grateful to DOE for the generous Grand Challenge Award 
which makes the calculations possible. 
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Figure Caption 

Fig. 1 The accumulated averages of the real part of the diagonal trace as estimated by 
the Z2 and Gaussian noises for three quark masses as functions of L, the number of noise 
vectors, are plotted in the left column. The right column shows the corresponding jackknife 
errors as functions of L. 

Fig. 2 The calculated axial form factor (7^(9^) and the proton electric form factor GE{q^) 
as a function of the momentum transfer —q^. The solid and the dash-dotted lines are fit to 
two experimental sets of data with different dipole masses. The dashed curve is the fit to 
the experimental Ge{q'^) with a dipole mass of 0.828 GeV. 

Fig. 3 The ttNN form factor g-KMNil'^) at the quark mass which corresponds to the physical 
pion mass. The solid/dashed line represents the monopole/dipole fit. 

Fig. 4 The ratio in eq. (|9|) for the scalar charge as a function of t for various quark masses. 
The slope is denoted by M and the P^r degree of freedom is also given. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9408007v3 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9408007v3 
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